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Abstract 

Recently, it was argued that new-physics (NP) effects in B decays can be approxi- 
mately parametrized in terms of a few quantities. As a result, CP violation in the 
B system allows one not only to detect the presence of new physics (NP), but also 
to measure its parameters. This will allow a partial identification of the NP, before 
its production at high-energy colliders. In this paper, we examine three methods 
for measuring NP parameters. The first uses a technique involving both b — ► s and 
b — > d penguin B decays. Depending on which pair of decays is used, the theoretical 
error is in the range 5-15%. The second involves a comparison of B — > itK and 
B — > 7T7T decays. Although the theoretical error is large (> 25%), the method can 
be performed now, with presently-available data. The third is via a time-dependent 
angular analysis of B — ► V\ V^ decays. In this case, there is no theoretical error, but 
the technique is experimentally challenging, and the method applies only to those 
NP models whose weak phase is universal to all NP operators. A reliable identifica- 
tion of the NP will involve the measurement of the NP parameters in many different 
ways, and with as many B decay modes as possible, so that it will be important to 
use all of these methods. 



1 Introduction 

Within the standard model (SM), CP violation is due to the presence of a complex 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix [1]. The principal goal 
of the study of CP violation in the B system is to test this explanation, and to 
find evidence for physics beyond the SM. As a result, much theoretical work has 
concentrated on signals of new physics (NP) in B decays [2]. 

At present, we have several experimental hints of new physics. First, within 
the SM, the CP asymmetry in B®(t) — ► J/ipK s should be approximately equal 
to that in decays dominated by the quark-level penguin transition b — > sqq (q = 
u,d,s). However, there is a 2.2a difference between the Belle measurement of the 
CP asymmetry in B®(t) — ► <pK s and that in B®(t) — ► J/ipK s [3]. In addition, BaBar 
sees a 3.0<j discrepancy between the CP asymmetries measured in B®(t) — ► i] f K s and 
Bj(t) — > J/ipK s [4]. Second, the latest data on B — » ixK decays (branching ratios 
and various CP asymmetries) appear to be inconsistent with the SM [5]. Third, 
within the SM, one expects no triple-product asymmetries in B — *■ <pK* [6], but 
BaBar has measured such an effect at 1.7<r [7]. 

It must be emphasized that none of these signals is statistically significant. Fur- 
thermore, the two experiments Belle and BaBar have not yet converged on any of the 
above measurements. Still, these signals are intriguing, particularly since the decays 
do share one thing in common: they all receive significant contributions from b — > s 
penguin amplitudes. If NP is indeed present, it is therefore plausible to suspect that 
it is the b — > s penguin which is principally affected. This is the assumption made 
in this paper. 

All new-physics effects in B decays are necessarily virtual. Thus, regardless of 
how evidence for physics beyond the SM is found — any of the above hints could give 
a statistically-significant signal of NP with more data — many models can explain 
any discrepancy. As a result, it has generally been assumed that the identification of 
the NP will have to wait until the new particles are produced directly at high-energy 
colliders. 

However, it was recently shown that this is not entirely true [8]. Briefly, the ar- 
gument is as follows. Following the experimental hints, we assume that new physics 
contributes significantly to those decays which have large b — > s penguin amplitudes. 
Consider now aB-»/ decay involving a b — > s penguin. The NP operators are 
assumed to be roughly the same size as the SM b — > s penguin operators, so the 
new effects are sizeable. At the quark level, the NP contributions take the form 
®np ~ sTibqTjq (q = u,d,s,c), where the Tjj represent Lorentz structures, and 
colour indices are suppressed. There are a total of 20 possible NP operators; each 
of them can in principle have a different weak phase. 

There are new-physics contributions to the decay B — > / through the matrix 
elements (/| O™ \B). Each of these can be written as 

{f\&™\B)=A k e<e< , (1) 



where <p q k and 8f. are the NP weak and strong phases associated with the individual 
matrix elements. However, it was argued in Ref. [8] that all NP strong phases 
are negligible compared to those of the SM. The point is that strong phases arise 
from rescattering. In the SM, this comes mainly from the tree diagram described 
at the quark level by b — ► sec 7 . However, note that this diagram is quite a bit 
larger than the b — > s penguin diagrams. That is, the strong phases associated with 
penguin amplitudes are due to rescattering from a diagram which is considerably 
bigger. However, the NP strong phases can come only from rescattering from the NP 
diagrams themselves, which are much smaller than the SM tree diagram. Thus, the 
generated NP strong phases are correspondingly smaller than their SM counterparts. 
That is, the NP strong phases are negligible compared to the SM strong phases. (A 
detailed discussion of small NP strong phases is presented in the Appendix.) 

Note that, in certain calculations of nonleptonic decays [9] it is claimed that the 
rescattering from the tree diagrams is negligible, but that annihilation terms, which 
are subleading (~ 0(1/to;,)), can be large. Large rescattering from the annihilation 
terms can generate a significant strong phase. If this scenario is true, then anni- 
hilation topologies associated with new-physics operators can also generate a large 
strong phase through rescattering. However, there is no general agreement on the 
importance of annihilation terms. Ultimately the size of annihilation diagrams is 
an experimental question, and can be tested by the measurement of decays such as 
B^ — ► D^D~ and B® — > K + K~ . If the annihilations terms turn out to be small, as 
expected from the 0(l/m&) suppression, then we can neglect the rescattering phase 
resulting from them. In our analysis we assume that annihilation-type topologies 
in the SM and with NP, which are power suppressed, are small, and therefore our 
argument that NP strong phases are negligible compared to those of the SM remains 
valid. 

The observation that the NP strong phases are negligible allows for a great sim- 
plification: one can now combine all NP matrix elements into a single NP amplitude, 
with a single weak phase: 

£</|0™|5>=.AV*«, (2) 

where q — u,d, s, c. Throughout the paper, we use the symbols A and $ to denote 
the NP amplitudes and weak phases, respectively. In the above, 

tan $ = £* A i Sin $ /gx 

E* A cos 4>1 ' 

Thus, all NP effects can be parametrized in terms of a small number of NP quantities. 
That is, we have an effective-lagrangian approach to new physics in CP violation in 
the B system. 



7 Thcre is also a f) ^ suu tree diagram, but it is described by the product of CKM matrix 
elements V* b V us , which is tiny. 



The argument the new-physics strong phases are negligible is quite general and 
applies to all NP models. However, this result can be obviated if special conditions 
are met. In particular, it does not hold if the NP is quite light, or if there is a 
significant enhancement of certain matrix elements. While these are disfavoured 
theoretically, the reader should be aware of these possible exceptions. 

Note that, in general, A q and Q q will be process-dependent. The NP phase $ g 
will be the same for all decays governed by the quark-level process b — > sqq only if all 
NP operators for the same quark-level process have the same weak phase. This is not 
uncommon. There are a number of NP models for which the weak phase is universal 
to all operators. These include models with Z- [10] and Z'-mediated [11] flavour- 
changing neutral currents (FCNC's), models in which the gluonic penguin operators 
have an enhanced chromomagnetic moment [12], and models with scalar- mediated 
FCNC's [13]. On the other hand, there are also NP models without universal weak 
phases, such as supersymmetric models with i?-parity breaking, left-right symmetric 
models and models with four generations. 

In Ref . [8] , it was shown that the A q and $ 9 can be measured using pairs of B 
decays which are related by flavour SU(3). One decay has a large b — > s penguin 
component and so receives new-physics contributions. The other has a b — »■ d 
penguin contribution. At present, there are no NP signals in processes which receive 
sizeable contributions from b — > d penguin amplitudes, such as B^ — > txtx . In this 
technique, and in others like it, we therefore assume that the NP does not affect 
decays involving b — ► d penguins. The measurements of the two decays permit the 
extraction of the NP parameters, which in turn allows one to discriminate among 
NP models and rule out many of them. We can thus partially identify the new 
physics, before high-energy colliders are used. 

In this paper, we provide a more detailed description of the method proposed in 
Ref. [8] to measure the NP parameters. We also examine two additional methods 
which can be used to obtain these NP parameters. 

The first new method involves B — > irK and B — ► ixix decays. Recently, it 
was shown that, within the SM, the full unitarity triangle can be extracted from 
measurements of B — »■ nK decays [14]. In order to do this, it is necessary to use 
flavour SU(3) to relate electroweak penguin operators to tree operators. On the other 
hand, if one assumes in addition the presence of new-physics b — > s amplitudes in 
B — ► iiK, it is straightforward to show that there is not sufficient information to 
extract the various SM and NP parameters. However, flavour SU(3) also relates 
B — > ttK to B — ► Till decays. Since the NP is not expected to affect these latter 
decays, one can use flavour SU(3) to obtain certain SM B — > ttK amplitudes from 
B — > 7T7T. With this information, it is possible to measure the NP parameters. The 
advantage of this method is that the analysis can be performed with present data; 
the disadvantage is that there is a theoretical error due to the assumption of flavour 
SU(3) symmetry. 

The second new method involves B — ► V1V2 decays, where V\ and V2 are vector 



mesons. These decays are very promising for finding evidence of physics beyond 
the SM. Suppose that the final state is such that (i) l / il / 2 = V1V2, and (ii) a 
single decay amplitude dominates in the SM. In this case, a time- dependent angular 
analysis of B(t) — ► V1V2 provides numerous signals of new physics [15]. Suppose 
further that a single NP amplitude is present, with a different weak phase from 
that of the SM amplitude and a (helicity-dependent) strong phase. The NP weak 
phase is assumed to be helicity-independent, which is the case for NP models whose 
weak phase is universal to all operators. In Ref. [15] it was shown that one can 
place lower bounds on the NP parameters. However, we have argued above that 
the NP strong phase is negligible, in which case the analysis is modified. As we will 
see, there are now more observables than theoretical parameters, so that one can 
measure the NP parameters in this system. Compared to Ref. [8], the advantage is 
that no theoretical input [flavour SU(3)] is required; the disadvantage is that it is 
difficult experimentally, and the analysis only holds for a certain class of NP models. 
We begin in Sec. 2 with a detailed discussion of the method involving b — > s 
and b — > d penguin decays which are related by SU(3). In Sec. 3 we turn to the 
analysis of B — »■ nK and B — ► tttt decays. The technique involving B — > V1V2 
decays is examined in Sec. 4. In Sec. 5 we discuss all three methods, stressing their 
advantages and disadvantages. All methods have their own unique features, and a 
complete analysis would ideally include all three techniques. We also examine two 
models of new physics, and show that different NP models lead to different patterns 
of NP parameters. This demonstrates that the measurement of the NP parameters 
does indeed discriminate among various models, and provides a partial identification 
of the NP. We conclude in Sec. 6. 



2 B Penguin Decays 

We begin with a description of the method proposed in Ref. [8] for measuring new- 
physics parameters. This technique closely resembles that of Ref. [16], which two 
of us (AD, DL) recently proposed for extracting CP phase information. Here the 
method is reversed. We assume that NP is present only in decays with large b — > s 
penguin amplitudes, and that the SM CP phase information is known (the SM 
phases can be measured using processes which do not involve large b — » s penguin 
amplitudes). We first study the general case in abstract terms. We then apply the 
method to specific decays. 

2.1 General Case 

We begin by considering a neutral B'° — > M[M 2 decay in which the final state M[M 2 
is accessible to both B° and B mesons. B'° can be either a B® or a B° s meson, and 
M[ and M 2 are two mesons. (If both M[ and M 2 are vector mesons, the final state 
can be considered as a single helicity state of M[M 2 .) This decay involves a b — ► s 



penguin contribution and is dominated by a single decay amplitude in the SM. (The 
case in which there are two significant SM amplitudes is discussed at the end of this 
subsection.) 

Since B'° — ► M[M 2 involves a6^s penguin amplitude, new physics is present. 
As discussed in the introduction, because the NP strong phases are negligible, the 
effect of the NP can be parametrized in terms of a single effective amplitude, with 
a NP weak phase. Thus, including the NP, the amplitude for B'° — ► M[M 2 can be 
written 

A(B'° -+M[M' 2 ) = A = A! ct e 1 ^ + A q J* q , (4) 

where A' ct and A q are the SM and NP amplitudes, respectively. Similarly, S' ct and 
$ g are the SM strong phase and NP weak phase, respectively. The amplitude for 
the CP-conjugate process, A, can be obtained from the above by changing the sign 
of the weak phase $ ? . 

The time-dependent measurement of B'°(t) — »■ M[M 2 allows one to obtain the 
three observables 

B ee 1 -(\A\ 2 + \A\ 2 ) = (A' ct ) 2 + (A q ) 2 + 2A ct A q cos5' ct cos%, 

a dir = i (\A\ 2 - \A\ 2 ) = 2A' ct A q sm8' ct sm$ q , (5) 

ai = Im (e" 2 *A*l) = ~{A ct ) 2 sin 2<p q M - 2A' ct A q cos 8' ct sin(20^ + $ 9 ) 

-(^) 2 sin(20t + 2$ 9 ). 

It is useful to define a fourth observable: 

a R = Re (e- 2j ^A*l) = (A ct f cos 20^ + 2A' ct A q cos ^ cos(20^ + $,) 

+ (^) 2 cos(20t + 2<l» g ). (6) 

The quantity a fl is not independent of the other three observables: 

a R 2 = B 2 - a dir 2 - a 2 . (7) 

Thus, one can obtain a R from measurements of B, a^r and o 7 , up to a sign ambiguity. 

In the above, <p q M (q = d, s) is the phase of B'°-B'° mixing. For B'° = B®, this 
phase is unaffected by new physics and thus takes its SM value, f3. The canonical 
way to measure this angle is via CP violation in B^(t) — > J/ipK s . However, there 
is a potential problem here: this decay receives NP contributions from O c NP ~ sbcc 
operators (the Lorentz and colour structures have been suppressed). On the other 
hand, the value of (3 extracted from B®(t) — > J/ipK s is in line with SM expectations. 
This strongly suggests that any O c NP contributions to this decay are quite small. 

The situation is somewhat different for B'° = B° s . In general, NP which af- 
fects b — > s transitions will also contribute to B®-B® mixing, i.e. one will have 
NP operators of the form sbsb. In this case, the phase of B° s -B® mixing may 



well differ from its SM value (~ 0) due to the presence of NP. The standard 
way to measure this mixing phase is through CP violation in B®(t) — > J/ipi] (or 
Bg(t) — > J J i/xf>). As with B%(t) — ► J/ipK s , these decays potentially receive 0^ p 
contributions. However, since the non-strange part of the i] wavefunction has a 
negligible contribution to {J/ipr]\O c NP \B®), this matrix element can be related by 
flavour SU(3) to (J/ipK s \ O c NP \B%) (up to a mixing angle). That is, both matrix 
elements are very small. In other words, we do not expect significant O c NP contribu- 
tions to B° s (t) — ► J/iprj, and the phase of B° s -B® mixing can be measured through 
CP violation in this decay, even in the presence of NP. 

Another decay which can be used to measure the phase of B®-B® mixing is 
B®(t) — > D+D~. Since the final state DfD~ is unrelated to J/ipt], it is logically 
possible that O c NP will have measurable effects in B®(t) — ► D+D~ . The easiest way 
to detect this is to measure the B®-B® mixing phase in both B®(t) — > J/ipi] and 
Bg(t) — > DfD~. If these two phases differ, this will clearly signal the presence of 
NP in b — > sec transitions. However, such measurements will not allow us to cleanly 
determine the magnitude and phase of O c NP . This can be done using the technique 
described in this section. 

Note that the expressions for a^ r and o 7 in Eq. (5) provide several clear signals 
of NP. Since B'° — ► M^M^ is dominated by a single decay in the SM, the direct 
CP asymmetry is predicted to vanish. Furthermore, in the SM the indirect CP 
asymmetry measures the mixing phase <p q M . Thus, if it is found that a^ r 7^ 0, or that 
4> q M differs from its SM value, this would be a smoking-gun signal of NP. Note also 
that, if it happens that the SM strong phases are small, a<n r may be unmeasurable. 
In this case, a better signal of new physics is the measurement of T-violating triple- 
product correlations in the corresponding vector- vector final states [6]. This is an 
example of the many NP signals present in B decays. However, these signals do not, 
by themselves, allow the measurement of the NP parameters. 

The three independent observables of Eqs. (5) and (6) depend on four unknown 
theoretical parameters: A q , A' ct , S' ct and $ g . Therefore one cannot obtain infor- 
mation about the new-physics parameters A q and $ g from these measurements. 
However, one can partially solve the equations to obtain 



(A' c 



2 = a R cos(20g / + 2$ g ) - q f sin(2^ + 2$ g ) - B 
cos2$ Q -l 



From this expression, we see that, if we knew A' ct , we could solve for $ g . 

In order to get A' ct we consider the partner process B° — > MiM 2 involving a 
b — > d penguin amplitude. In the SM this decay is related by SU(3) symmetry to 
B'° — > M^M^. (In some cases, this relation only holds if one neglects annihilation- 
or exchange-type diagrams [16], which are expected to be small.) B° can be either 
a B" or a B® meson and, as with B'° — > M^M'^ it is assumed that both B° and 
B° can decay to the final state MiM 2 - The partner process can be a pure penguin 
decay, or can involve both tree and (non-negligible) penguin contributions. 



Since b — > s transitions are not involved, the amplitude for B° — > M\M 2 receives 
only SM contributions, and is given by 

A(B°^M 1 M 2 ) = A u v: b V ud + A c V; b V cd + A t V t * b V td 

= (A u -A t )V: b V ud + (A c -A t )V; b V cd 

= Aute^e^+Acte^ , (9) 

where A ut = \(A U - A t )V* b V ud \, Act = I (A c - A t )V* b V cd \, and we have explicitly 
written the strong phases 5 ut and S ct , as well as the weak phase 7. In the above, we 
adopt the c-quark convention [17], in which CKM unitarity is used to eliminate the 
t-quark term. 

As with B'° — >■ M^M^, the time-dependent measurement of B°(t) — ► MiM 2 
allows one to obtain three independent observables [Eqs. (5) and (6)]. These ob- 
servables depend on five theoretical quantities: Act, A u t, $ = but ~ b ct) 7 and the 
mixing phase (p q M . However, as discussed above, (p q M can be measured independently 
using processes which are unaffected by new physics in b — > s penguin amplitudes. 
The weak phase 7 can be measured similarly. For example, it can be obtained from 
B^ — ► DK decays [18]. Alternatively, the angle a can be extracted from B — ► tttt 
[19], B — > pn [20] or B — > pp decays [21], and 7 can be obtained using 7 = 7r — (3 — 7. 
Given that these CP phases can be measured independently, the three observables 
of B°(t) — > MiM 2 now depend on three unknown theoretical parameters, so that 
the system of equations can be solved. 

In particular, one can obtain Act' 

2 = a R cos(20^ f + 27) - q I sin(20 M + 27) - B 
cos 27 — 1 

where a R , a 7 and B are the observables found in B°(t) — >■ M1M2. 
The key point is that, in the SU(3) limit, one has 

Act = \A' ct , (11) 

where A = 0.22 is the Cabibbo angle. With this relation, the extraction of Act 
from B°(t) — > MiM 2 yields A' ct . Thus, Eq. (8) can be used to solve for the new 
physics phase $ 9 . The NP amplitude A q can also be obtained. There is a theoretical 
error in the above relation due to SU(3)-breaking effects. However, various methods 
were discussed in Ref. [16] to reduce this SU(3) breaking. All of these methods 
are applicable here. In the end, depending on which pair of processes is used, the 
theoretical error can be reduced to the level of 5-15%. 

At this point we can make an important general observation. As noted in the 
introduction [Eq. (3)], the new-physics weak phase $ g depends on the matrix ele- 
ments of the various NP operators for the particular process considered (Ai), as well 
as the corresponding weak phases </>f. Thus, in general, the value of $ g extracted 
from two distinct decay pairs with the same underlying b — > sqq transition will be 



different. There are two reasons for this. First, certain operators which contribute 
to one process may not contribute in the same form to another. (For example, one 
decay might be colour-suppressed, while the other is colour-allowed.) Second, even 
if the same operators are involved (with the same form) in two b — > sqq decays, the 
matrix elements of the various operators will depend on the final states considered. 
Thus, the A* in Eq. (3) are process-dependent in general, and the value of the phase 
$ g depends on the particular decay pair used. However, if all NP operators for the 
quark-level process b — > sqq have the same weak phase <p q , then the NP phase Q q 
will be the same for all decays governed by the same quark-level process. Hence it 
is important to measure the phase Q q in more than one pair of processes with the 
same underlying quark transition. If the effective phases are different then it would 
be a clear signal of more than one NP amplitude, with different weak phases, in 
b — »■ sqq. Furthermore, in some NP models, the phases for the different underlying 
quark transitions b — > sqq are related, so that the NP phase is independent of the 
quark flavour. 

In the above method, we have assumed that the decay B'° — > M' X M' 2 is dominated 
by a single decay amplitude in the SM. This is the case only for the quark-level 
decays b — > sqq (q = d, s, c). However, it is straightforward to adapt this technique to 
b — > suu, for which B'° — > M' X M ! 2 receives both tree and b — > s penguin contributions 
in the SM. The process B® — > K + K~ is an example of such a decay. Including the 
new-physics contribution, the amplitude for such decays can be written 

A(B'° -► M[M' 2 ) = A' ut e i7 e w «« + A ct e iS « + A u e 1 ^ . (12) 

Here, assuming that 7 and the mixing phase <p q M are known, the three independent 
observables in this decay depend on six unknown parameters: A' ct , A' ut , S' ut , 5' ct , A u 
and $ u . In this case, in order to solve for the NP parameters, one needs three pieces 
of information. These can be obtained as follows. Measurements of the partner 
process allow one to extract Act, A u t an d 5 = S u t — S c t. We now assume that 

A ct = XA' ct , \A ut = A' ut , S' = 8, (13) 

where 5' = S' ut — S' ct . These assumptions then permit the extraction of A u and $ u 
from measurements of B'° — > M' X M' 2 . 

However, the theoretical uncertainty here due to SU(3) breaking is considerably 
larger than in the case where B'° — > M^M'^ is dominated by a single amplitude in 
the SM. Not only do we relate two amplitudes instead of one [Eq. (11)], but we also 
assume that two strong phases are equal. Thus, the NP parameters A u and $ u can 
be obtained in this way, but we expect a larger theoretical error. 

2.2 Specific Decays 

In Ref. [16], we showed that there are twelve decay pairs B° — > M1M2 and B' — >■ 
M[M2 which can be used to obtain CP phase information in the SM, with a small 
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NP Parameters 


B'°(t) - 


^M[M' 2 


B°(t) - 


-> MiM 2 


$c A c 


B° s (t) - 


- DtD- 


B° d (t) - 


^ D+D- 


® S ,A S 


B° d (t) ~ 
B° s (t) - 




B° s (t) - 

b°M - 




® d ,A d 


B° s (t) - 
B° s (t) ~ 
B° d (t) - 
B° d (t) - 


-> K*°p° 


B° d (t) - 
B° d (t) - 
B° d (t) - 

b°M - 


^p°p° 
■+ K*°p° 


®u,A u 


B° s (t) - 


^K+K- 


B° d (t) - 


->■ 7T + 7r _ 



Table 1: For each set of new-physics parameters ^4 9 and $ g , we list the 5 decays 
(B'°(t) — >• M{Mg) and their partner processes (B°(t) — ► MxM 2 ) which can be used 
to measure them. 



theoretical error. (In fact, there are more, since many of the particles in the final 
states can be observed as either pseudoscalar (P) or vector (V) mesons.) Many of 
these decay pairs can also be used to measure the new-physics parameters, assuming 
that the NP contributes significantly only to the b — » s decays, and that the SM CP 
phases have already been measured using non-6 — > s processes. 

As noted earlier, assuming that new-physics strong rescattering is negligible rel- 
ative to that of the SM, all NP effects can be parametrized in terms of the effective 
NP amplitudes A q and weak phases $ g (q = u,d,s,c), independent of the type of 
underlying NP [8] . Of the above twelve decay pairs, seven involve only neutral B- 
mesons and have final states M' X M' 2 accessible to both B° and B . These can be used 
to measure the NP parameters A q and $ g (q = d, s, c) using the method outlined in 
Sec. 2.1. For A u and $ u , we have to use a decay pair in which B' — ► M' X M' 2 receives 
both tree and penguin contributions in the SM. There is one such possibility. Thus, 
all four sets of NP parameters can be obtained from measurements of pairs of B 
decays. 

The decay pairs are listed in Table 1, along with the new-physics parameters 
probed. We have several comments about these. 

There are three reasons why certain decays are written in terms of vector- vector 
(VV) final states, while others involve pseudoscalar-pseudoscalar (PP) states. First, 
some decays involve a final-state 7r°. However, experimentally it will be necessary 
to find the decay vertices of the final particles. This is virtually impossible for a 
7T°, and so we always use a p° [22]. Second, some pairs of decays are related by 
SU(3) in the SM only if an (ss) quark pair is used. Unfortunately, there are no 
P's which are pure (ss). The mesons i] and rj' have an (ss) component, but they 
also have significant (uu) and (dd) pieces. As a result the decays B' — > M[M 2 and 
B° — > M\M 2 are not really related by SU(3) in the SM if the final state involves 



an r) or r/. We therefore consider instead the vector meson <fi which is essentially a 
pure (ss) quark state. Finally, we require that both B° and B° be able to decay to 
the final state. This cannot happen if the final state contains a single K° (or K°) 
meson. However, it can occur if this final-state particle is an excited neutral kaon. 
In this case one decay involves K*°, while the other has K*°. Assuming that the 
vector meson is detected via its decay to K s 7r° (as in the measurement of sin 2/3 via 
Bj(t) — ► J/ipK*), then both B° and B° can decay to the same final state. 

Apart from these three restrictions, the final-state particles can be taken to be 
either pseudoscalar or vector. Indeed, it will be useful to measure the NP parameters 
in modes with PP, PV and VV final-state particles, since different NP operators are 
probed in these decays. For example, within factorization, certain scalar operators 
cannot contribute to PV or VV states if their amplitudes involve the matrix element 
(V\ q^L^Q |0). In general, the matrix element of a given operator will be different 
for the various PP, PV and VV final states. Thus, the measurement of the NP 
parameters in different modes will provide some clues as to which NP operators are 
present. 

In addition, if it is found that $ 9 is different for decays governed by the same 
underlying quark-level transition, it will indicate the presence of more than one NP 
amplitude, with different weak phases. 

For the NP parameters A q and $ 9 (q = d, s, c) , the theoretical error is due to 
SU(3) breaking in Eq. (11), and can be reduced to the range 5-15% [16]. On the 
other hand, the only way to measure $ M and A u is to use B®(t) — *■ K + K~ and 
Bj(t) — ► 7r + 7r~. However, B® — > K + K~ has both tree and penguin contributions. 
In order to obtain $ M and A u , it is therefore necessary to make the three assumptions 
in Eq. (13). In the context of measuring the angle 7, the SU(3) breaking in B® — ► 
K + K~ and B® — ► 7r + 7r~ was examined in Refs. [23, 24]. In the framework of naive 
factorization or QCD factorization [25], it can be shown that, as long as annihilation- 
type topologies are small, the double ratio of amplitudes \ 2 (A' ct /A' ut )/(A c t/A u t) has 
small SU(3) breaking. However, Eq. (13) does not involve this double ratio of 
amplitudes - it involves single amplitude ratios (and an equality of strong phases). 
In this case, even within QCD factorization with small annihilation-type topologies, 
there are several sources of SU(3) breaking which are not under total control. The 
SU(3) breaking comes from the difference between unknown B° s —> K and B® — > 71 
form factors, differences in the light cone distributions of the kaon and the pion, 
and other subleading but potentially important unknown soft physics [26]. As a 
result, putting all these SU(3)-breaking effects together, the theoretical error in the 
extraction of $ M and A u is quite a bit larger than for the measurement of the other 
NP parameters. 

Note that only one pair in Table 1 involves only B^ decays. The others will 
require the time-dependent measurement of B° s decays. However, this will be difficult 
experimentally, as B®-B® mixing is large. For this reason the decay pair B®(t) — > 
K*°p° and B^(t) — >■ p°p° may be the most promising for measuring NP parameters 
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using this method. 

3 B — > 7rK and B — ► mr Decays 

In this section we consider B — > 7iK and B — ► tttt decays. It is well known that 
it is possible to express the amplitudes for B decays to two pseudoscalars in terms 
of a number of distinct SU(3) operators. This is equivalent to a description in 
terms of diagrams [27]. Neglecting the exchange- and annihilation-type diagrams, 
which are expected to be small for dynamical reasons, but including electroweak 
penguin contributions (EWP's), there are five diagrams [28]: (1) a colour-favored 
tree amplitude T (or T"), (2) a colour-suppressed tree amplitude C (or C), (3) 
a gluonic penguin amplitude P (or P'), (4) a colour- favored electroweak penguin 
amplitude P EW (or P' EW ), and (5) a colour-suppressed electroweak penguin amplitude 
P° w (or P' EW ). In the following, we denote all diagrams contributing to b — > d (b — *■ s) 
decays without (with) primes. 

As described in Sec. 2, the penguin diagram actually contains several pieces: 

p = p u v: b v ud + p c v; b v cd + p t v t * b v td 
= (p u -p t )v: b v ud + (p c -p t )v; b v cd 

= P ut e l ^ + P ct . (14) 

In the above, we have explicitly written the weak phase 7; the amplitudes implicitly 
include the strong phases. As in Sec. 2, we have adopted the c-quark convention 
[17] in passing from the first line to the second. For b — ■> s decays, we can write P' 
analogously to the above, except that we expect |P^I ^ \P'ct\ srnce iKTb^ws/Kb^csl — 
2%. In Sec. 2 we neglected the P' ut term. We will eventually do something similar 
here as well, but for the moment we keep all terms in the B — > irK amplitudes. 

For b — > d decays, the EWP contributions are expected to be negligible. How- 
ever, they are important for b — ► s transitions [28]. It was recently shown that, 
to a good approximation, the EWP's can be related to tree operators using Fierz 
transformations and SU(3) symmetry [29]. Ignoring exchange- and annihilation-type 
diagrams once again, the relations are 

P' EW = \^±^R{T' + C') + \^^R{T'-C), 
4 ci + c 2 4 ci - c 2 

P' E C W = 3 - C -^-^R(T' + C)- 3 - C -^^R(r-C'), (15) 

4 ci + c 2 4 ci - c 2 

where the q are Wilson coefficients. Here, the weak phases have been factored out, 
so these relations include only strong phases. In the above, 



R 



v; b v ts 



v: h v us 
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(16) 



where p and r) are the CKM parameters. (Note that the three CP phases in the 
unitarity triangle are functions of p and 77.) Thus, all B — > PP amplitudes can 
be expressed in terms of T, C, P u t, P c t (or their b — > s equivalents) and the weak 
phases. 

We now turn to the four B -»■ nK decays B + -»■ vr+fC , 5+ -»• 7r°^+, 5^ -> 
7r^7^ + and B d — > tt ^ . There are 9 measurements that can be made of this sys- 
tem: four branching ratios, four direct asymmetries, and one indirect asymmetry 
(in Bj(t) — ► 7i° K s ). However, assuming that the P' ut term is negligible, in the SM 
the amplitudes can be expressed in terms of 7 parameters: three diagram magni- 
tudes, two relative strong phases, and two CKM parameters. Thus there is enough 
information in B — > nK decays to reconstruct the full unitarity triangle [14] . 

We now consider new physics. B — > nK decays are b — > s transitions, so that 
NP can affect these decays. There are two classes of NP operators, differing in their 
colour structure: s a Tib a qpTjqp and s a Tib/3q/3Tjq a . The first class of NP operators 
contributes with no colour suppression to final states containing qq mesons. The 
second type of operator can also contribute via Fierz transformations, but there is 
a suppression factor of 1/N C , as well as additional operators involving colour octet 
currents. Similarly, for final states with sq mesons, the roles of the two classes of 
operators are reversed. We denote by A' ,q e 1 ®* and A' c ' q e l ®i the sum of NP operators 
which contribute to final states involving qq and sq mesons, respectively. Here, $' 
and $' c are the NP weak phases; the strong phases are zero. We stress that, despite 

the "colour-suppressed" index C, the operators A rc ' q e « are not necessarily smaller 
than the .A'-V*'. 

Including these NP operators, the B — > nK amplitudes can be written 

A(B+^n+K°) = P' ut e^ + P' ct - l -P' E c w + A^ d e^ C , 
V2A(B + -► 7T° K + ) = -P' ut e 17 - P' ct - T' e* 7 - C e* 7 - P' EW 

_ lp'c w _ A'' u e^ + A'' d e^ - A' c ' u e l ^ C , 

A(B° d - n-K+) = -P' ut e* 7 - P' t - T> e* 7 - h>'£ w - A' c ' u e^° , 

V2A(B° d -> n°K°) = P' ut e^ + P' ct -C'e^-P' EW - l -P' E c w 

- A'' u e 1 ^ + AVi + A ,c4 e^ . (17) 

Here, each of P' uV P' ct1 T' and C include a (different) strong phase. Note that 
A'' u e u and A' ,d e l ®d always appear in the same combination above. We therefore 
define A! ,comb e % ^ > = —j\>> u e l K _|_ A'^e 1 ®*. It is not possible to distinguish the two 
component amplitudes. 

In the presence of NP, the amplitudes can be written in terms of 16 theoretical 
quantities: 7 amplitude magnitudes (\P' ut \, \P' ct \, \T'\, \C'\, \A'> comb \, \A' c ' d \ and 
|^4 /c '"|), 4 relative strong phases, 2 SM weak phases, and 3 NP weak phases. (In 
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the following, we generically refer to the strong phases, SM weak phases and NP 
weak phases as "<5", "0" and "<&," respectively.) Since we have 16 parameters and 
only 9 measurements, it is clear that we cannot measure the NP parameters using 
B — ► ttK alone. It does not help to make the approximation that P' ut is negligible. 
Important information can be obtained from measurements of the B — ► tttt 
system. As per our assumptions, new physics does not affect such decays. Neglecting 
EWP contributions, which are expected to be small, the B — > tttt amplitudes can 
be written 

V2A(B + -► tt+tt ) = -T e i7 - C e %1 , 

A(B° d - 7T+7T-) = -Te^-P ut e^-P ct , 



v^A(SS -^ ttV) = -CV 7 + P ut e i7 + P 



ct 



The indirect asymmetry in P°(i) — > 7r + 7r _ also involves the phase of B^-B^ mixing, 
(3, so the 6 B — »■ 7r7r measurements [three branching ratios, two direct asymmetries, 
and one indirect asymmetry (in B®(t) — »■ 7r + 7r~)] are a function of 7 theoretical 
parameters. However, if one assumes that /5 is measured in B®(t) — »■ J/ipK s , all 
theoretical quantities can be extracted [30]. 

Within flavour SU(3) symmetry, the P — ► 7r7r amplitudes are related to those in 
B^ttK: 

T _ C _ P ut _ X Pd _ P EW _ P EW _ X . . 

T C P ut 1 - A 2 /2 ' P c ' t P^ w P>c w 1 - A 2 /2 ' l ) 

where the various amplitudes include the strong phases. With these relations, we 
can combine the information obtained in B — > tttt and B — >■ 7r7^ decays. As detailed 
above, there are a total of 16 theoretical parameters. However, there are now also 
16 experimental measurements: 6 in B — > tttt, 9 in B — > ttK and the extraction of 
/3 in B®(t) — > J/ipK s . It is therefore in principle possible to solve for all theoretical 
unknowns, and we would thus measure the NP parameters. 

Unfortunately, solving 16 nonlinear equations in 16 unknowns will lead to a 
large number of discretely-ambiguous solutions. When one adds the experimental 
errors, the solutions will be smeared out, and the values of the NP parameters will 
essentially remain unknown. To remedy this, we adopt the procedure of Ref. [31] in 
order to reduce the number of theoretical unknowns. First, we neglect the P' ut term 
in the decays B + — > 7r + fT and B + — > tt°K + . Second, we remove the dependence 
on P ut by redefining T and C: 

f = T + P ut , C = C-P ut , (20) 

with similar redefinitions for the primed quantities. Finally, the relations in Eq. (15) 
no longer hold when T' and C are used. We therefore neglect the amplitude P'£ w . 
With this, there is a single relation between P' EW and the tree diagrams: 

P> EW = 3 - C ^±^R(f + C>). (21) 

2 ci + c 2 
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Since \P^ t \, \P'e W \ ^ \P'ct\ we ex P ec t the error associated with these approximations 
to be small. Indeed, in Ref. [31], SM fits both with and without the approximations 
were performed, and little difference was found. 

With the above approximations, the B — ► itK amplitudes take the form 



A(B+ - 
V2A(B+ - 

V2A(B° d 



■ 7T°K + ) 



P' 



A 



IC.diQ 



IC 



-PL - T' e i7 - C e i7 - P' + A'' comb e^ - A 



\IC,uJ,&,9 



ct 

-P' 

r ct 



T' e* 7 - A c ' u e^ 



P' ct -C'e^-P' EW + A 



l,comb i& 



+ A 



IC,d i$ 



(22) 



while those for B 



7T7T are 

V2A{B + 
MB° d 

V2A(B° d -. 7T 7T ) 



7T + 7T ) 
0_0\ 



-T e 47 - C e 17 

-fe l7 -P ct , 
-Ce^ + P ct . 



The amplitudes with tildes are related as in Eq. (19): 






C 



A 



1 - A 2 /2 ' 



(23) 



(24) 



There are now 14 theoretical unknown quantities, but 16 measurements. Thus, we 
can solve for the NP parameters with few discrete ambiguities. (In practice, one will 
fit for all parameters.) 

Of course, we have assumed perfect SU(3) symmetry in this procedure [Eqs. (19) 
and (24)]. However, we know that there may be significant SU(3)-breaking effects. 
In Ref. [31], it was noted that factorization appears to hold for colour-allowed tree 
diagrams, so that \T'/T\ ~ f K / f w , but that the SU(3) breaking in the other relations 
would be left to experiment. In our case, NP is present in all b — > s decays, and this 
will mask any SU(3)-breaking effects. One possibility is to assume that the ratios of 
magnitudes of amplitudes are known, but no assumption is made about the strong 
phases. That is, we write 



V 



/: 



T 



~c 



/« 



c 



P 



ct 



P' 



I 



p ■ 



(25) 



The quantities fr, fc an d fp are calculated using some theoretical model (e.g. 
QCD factorization [25]), but all strong phases are taken to be additional theoret- 
ical unknowns. The problem here is that this adds two theoretical quantities to 
the procedure (two strong phases in B — > im decays), so we once again have 16 
measurements and 16 unknowns. As discussed above, this leads to a large number 
of discretely-ambiguous solutions. For this reason, it is probably best to assume 
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that the strong phases of primed and unprimed amplitudes are equal, as with per- 
fect SU(3) symmetry, and that the magnitude ratios are given by Eq. (25). In this 
case the above procedure will yield the NP parameters, but with sizeable theoretical 
errors. 

Above, we have concentrated on B — ► PP decays, where P is a pseudoscalar. 
However, the analysis holds equally for B — > VV decays (V is a vector meson). In 
this case, an angular analysis must be performed. Note that we have argued that 
exchange- and annihilation-type contributions to the B — ► PP decays are expected 
to be negligible. However, in some approaches to hadronic B decays, such amplitudes 
may be chirally enhanced if there are pseudoscalars in the final state [9, 25]. On 
the other hand, such chiral enhancements are not present for VV final states, so 
this is a potential point in favour of B — > VV decays. It has been recently claimed 
that annihilation terms can also be big in certain B — > VV decays in spite of the 
lack of chiral enhancement [32]. As noted earlier, ultimately, the size of exchange 
and annihilation diagrams is an experimental question, and can be tested by the 
measurement of decays such as B d — > DfD~ and B® — > K + K~ . 

One can even apply the method to PV final states, but things are considerably 
more complicated in this case since the NP contributes differently to the PV and 
VP final states. In this case, one must correspondingly increase the number of NP 
parameters: there are now a total of 30 parameters in e.g. B — * ttK* and B — > pK 
decays. These parameters are to be fitted to 32 measurements (9 in B — *■ ttK*, 9 
in B — > pK, 13 in B — ► p7r, (3 from B d (t) — > J/ipK s ). Thus, while one can solve for 
the NP parameters in principle, in practice the analysis is very complicated. 

We can do better with the B — > nK/B — > tttt method if we perform a semi- 
mo del- independent analysis by making an assumption about the general form of 
A n s. We illustrate this below. 

3.1 Isospin-conserving new physics 

One possibility is to make the general assumption that the new physics is isospin- 
conserving. For example, this occurs in NP models in which the gluonic penguin 
operators have an enhanced chromomagnetic moment [12]. In this case, 

ji^^'v, — A'^^K ? A' c ' u e i ^ = A' c ' d e i<!> '° = ,4"V*' C . (26) 

This in turn implies that A!' comh e 1 ^ is zero, so that there is just one A' remaining. 
The B — > ttK amplitudes now take the form: 

A(B + ^n + K°) = P' ct + A' c e^' C , 

V2A(B + -► tt° K+) = -P' ct - f e* 7 - & e* 7 - P' EW - A /c e^' C , 

A{B° d -► 7i- K+) = -P' ct - f e i7 - A c e l *' C , 

V2A(B° d -. n°K°) = P' ct - C e^ - P' EW + A lc e w ° . (27) 
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There are now only 10 theoretical parameters: 4 amplitude magnitudes, 3 <5's, 2 0's 
and one $. Recall that there are 9 measurements in B — » 7T.K" alone. We therefore 
need only one additional measurement to be able to extract all parameters, including 
those related to NP: A' c and $ /c . Ideally, in order to reduce discrete ambiguities, 
we would have two additional measurements. These can come from independent 
measurements of the SM phases (the 2 0's). Thus, for this particular type of new 
physics, we do not need measurements in the B — > ttti system at all. 

Indeed, it is not necssary to make any assumptions about the absence of new 
physics in decays with b — *■ d penguins, though we must assume that the phase of 
B®-B® mixing is measured in B d (t) — ■> J/ipK s . If we assume that NP is also present 
in B — > 7T7T decays, these amplitudes take the form 



f t r~< A otfyC 



C,d i® 



\PlA(B + -► tt+tt ) = -f e i7 -Ce n + A comb e^ - A c ' u e^ + A 

A(B° d -► tt+tt-) = -f e^ - P ct - A C '"J*° , 
V2A(B° d -► 7T°7r ) = -C e' 7 + P ct + A comb e^ + A c ' d e^ , (28) 

where 

A comb e i<5> = _ju e i$ u + ^d e i$ d _ ^ 

If the NP is isospin-conserving, this implies that 

,AV*« = .AV*" , A C 'V*« = ,A C >V*? = AV*"^ . (30) 

In this case, A comb e r ^ vanishes, and we are left with only one A, namely A c e 1 ^ ' . 
Assuming perfect SU(3) symmetry, the B — ► nK and 5 — ► tttt amplitudes are 
described by a total of 12 theoretical quantities. With 16 measurements, it is possible 
to extract all parameters, including those describing the NP in the b — » s and b — » d 
transitions. 



4 B ^V{V 2 Decays 

In this section we examine Z? — ► V1V2 decays in which l / il / 2 = ^1^2- We consider 
decays which are described by the quark-level transitions 6 — ► ccs, b — ► sss, or 
6 — ► 5<i(i. Within the SM such decays are dominated by a single weak decay ampli- 
tude, and their weak phase is essentially zero in the standard parametrization [1]. 
(As noted above, there are two significant contributions - - the tree and penguin 
amplitudes -- for decays described by b — ► suu.) Since these are all b — > s transi- 
tions, there are new-physics contributions. Note that, in this method, we make no 
assumptions about NP in b — > d decays. 

Suppose that the underlying new-physics model is such that the weak phase is 
universal to all NP operators. As discussed in the introduction, this holds for a 
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large number of NP models. In this case, the NP weak phase Q q will be helicity- 
independent. Taking into account the fact that the NP strong phase is negligible, 
the decay amplitude for each of the three possible helicity states may be written as 



A x = Amp(B - 


+ W 2 )a = 


= axe^+A\e^ 


A x = Amp(B - 


^V 1 V 2 ) X = 


= axe** + Ale-**- 



(31) 

where ax and A\ represent the helicity-dependent SM and NP amplitudes, respec- 
tively, the 5x are the SM strong phases, and the helicity index A takes the values 
{0, ||, _L}. Using CPT invariance, the full decay amplitudes can be written as 

A = Amp(B -► V X V 2 ) = A g + A im + i A ±g± , 

A = Amp(B -► VM) = A g + A m - i A ±g± , (32) 

where the gx are the coefficients of the helicity amplitudes written in the linear 

polarization basis. The gx depend only on the angles describing the kinematics [33]. 

Using the above equations, we can write the time-dependent decay rates as 

r$(t) - V 1 V 2 ) = e- rt J2 (Aa. ± S Act cos(AMt) T Px„ sm(AMtj)g x g ff . (33) 

A<cr 

Thus, by performing a time-dependent angular analysis of the decay B(t) — ► ViV 2 , 
one can measure 18 observables. These are: 

^xx = \(\Ax\ 2 -\Ax\ 2 ), 

A|| = Re(A||A*+A||A*), 

Z ll0 = Re(A ll A*-A ll A*), 

p ±± = Im(e- i4> M A* ± Aj 

Plio = -Im(e-^M [AlA + A* A^}) , Pii = -Im^e^M A* A t ) , (34) 

where % = {0, ||}. As in Sec. 2.1, <p q M is the weak phase factor associated with Bq-B® 
mixing. Note that the signs of the various pxx terms depend on the CP-parity of the 
various helicity states. We have chosen the sign of p„ to be — 1, which corresponds 
to the final state <pK*. 

For measuring new-physics parameters, the key point is the following. There are 
a total of six amplitudes describing B — ► V\V 2 and B — > V\V 2 decays [Eq. (31)]. At 
best one can measure the magnitudes and relative phases of these six amplitudes. 
Thus, of the 18 observables, only 11 are independent. However, these observables 
are a function of only 11 theoretical parameters 8 : three a A 's, three *4. A 's, (f> q M , $ g , 

8 If the NP weak phase is assumed to helicity-dependent, then there are more theoretical pa- 
rameters than there are measurements, and we cannot solve the system of equations. 
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A A a = 


^(I^a| 2 + MU| 2 ), 


A_u = 


-Lm(A ± A*-A ± A* 


£_u — 


-Im(A ± A*+A ± A* 


= Re(e" 


-^[AlAi+A^Aj 


— Imfe 


- i( t>M[A\A Q +AlA\\\ 



and the three strong phases 5\. In addition, as discussed in Sec. 2.1, (p q M can be 
measured independently, so we effectively have 11 equations in 10 unknowns. The 
solution will have discrete ambiguities, but many of these can be removed using the 
additional observables. 

Thus, as advertised, if new physics is found, it is possible to measure the NP 
parameters via a time- dependent angular analysis of B — ► V1V2 decays. To be 
specific, A s x and $ s can be extracted from B^ — ► <f)K*° or B® — > 00. The decays 



B® — ► i^*°p° and 5° — ► K*°K*° can be used to measure ^ and $ d . (In the 
> d —r K*°p°, there is a small theoretical error due to the neglect of the 



decay B® 

colour-suppressed b — > sufi tree contribution.) Finally, measurements of the decay 

5° — > D* + D*~ can be used to obtain A c x and $ c . 

Note that this analysis is done within the context of a single b — > s -B — »■ V1V2 
decay. In this case, as in the Z? — > ttK/B — ► 7r7r method with isospin-conserving 
new physics, no assumption is necessary about the absence of new physics in decays 
with 6 —> d penguins. The only assumption needed is that the phase of B®-B® 
mixing, which may be affected by NP, can be extracted from B®(t) — * J/ipK s or 
B° s (t) -> J/^. 



4.1 Explicit Solution 

Under the assumption that 0^ is known independently, we can construct an analytic 
solution (this follows closely the analysis of Ref. [15]). In terms of the theoretical 
parameters, the explicit expressions for the observables are as follows: 



Aaa 
£aa 

A l|0 

Silo 



a\ + {A q x f + 2a x Al cos ^a cos $ g , 

2axA q x sin 5a sin $ g , 

2 [a_i_v4f cos 5^ — aj^4^ cos 5,] sin $ g , 

2 oyoo cos(<5|| — So) + a||^4.Q cos ^|| cos &q + a oA q , cos 5o cos $ q 
—2 [a^aj sin(5j_ — 5j) + a±Aj sin 5j_ cos $ g — a^^ sin <5j cos $ 



,4iU2 



?J • 



0||^4q s i n ^11 + a o-A q sin So 
1 



sin$ 



q ■ 



Pii = a\ sin 20^ + 2a t A q cos ^ sin(20« f + $,) + (^) 2 sin(20« f + 2$,), 
P±± = -ai sin 20^ - 2a ± ^l cos S± sin(20« f + $,) - (A q ± ) 2 sin(20^ + 2$ 
P±i = 2 [aja^ cos(5j — S±) cos 20* f + a±A q cos 5^ cos(20^ + $ g ) 
+ ai^4i cos £ cos(20« f + $,) + A ? ^I cos(20^ + 2$,)] 
2 



?/• 



Plio 



OoO|| cos(<5|| — <5 ) sin 20^ + a||^4.Q cos <5|| sin(20^ + <3> g ) 
+ oo^H cos 5 sin(20« f + $,) + A%A\ sin(20« f + 2$ g ) 



(35) 
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For B — > V{V<2 decays, the analogue of the usual direct CP asymmetry a^-f is the 
helicity-dependent quantity a x r = Eaa/Aaa- We define the related quantity 



y x = y/l - S| A /AL • (36) 

Similarly, the value of sin 20^ measured in B — ► V1V2 decays can depend on the 
helicity of the final state: p\\ can be recast in terms of a measured weak phase 
(2(/>^)™ eas , defined as 

sin 2 (^)r = /A2 ±PA " , (37) 

V A AA - S AA 

where the + (— ) sign corresponds to A = 0, || (_L). 

Using the expressions for Aaa, £aa and ((f) q M )^ eas above, one can express a x and 
A\ as follows [15]: 

2a^sin 2 $ g = A AA (1 - yx cos(9 x - 2$,)) , 
2(^) 2 sin 2 $, = A A a(1-2/acos#a) , (38) 

where 9 X = (2(f) q M ) x neas — 2(p q M . Using these expressions, along with those for Aaa and 
Saa [Eq. (35)], we can solve for tan<^: 

tan5 = S AA sin 2$, 

^ A A X x(-2cos 2 <5> q + y x cos(9 x -2<5> q )+y x cos9 x )' { ' 

This equation expresses 5 X in terms of observables and <& q . 

How one proceeds further depends on which other observables are available. 
Suppose that Aj_j and Sj_j have been measured. These two observables can be 
expressed as [15] 

£_u = PiP±[(£±cn-€iCr±) cosAi-(l + £i£ ± + <7i a^) sin Ai] , 

A_u = PP± [(£l - &) cos A, - (a t + a ± ) sin A, ] , (40) 



where A,- = S 1 — Si, and 



i; 



P 2 = Aaa(1 - yx cos(20 A - 2%)) , 

t = Aaa1/a sin(2fl A -2$ g ) 

SA — p2 ) 

^A 

-A ES %. (41) 



^ 



Eqs. (40) can be solved for A3 



tan Ai = n 1 c c J, r* v — ? — 7 V • 42 ) 

(1 + Ki £1 + <Ti crj_) Aj_j - E ±i (<jj + ct_l) 
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This expresses tan A, in terms of observables and $ ? . However, we can also write 

tan 5± — tan & . , „. 

tanA t = / * . (43) 

1 + tano^ tanoj 

Eqs. (39), (42) and (43) can then be combined to give a single equation as a function 
of $ g . This can be solved to get the new-physics weak phase, which will permit the 
measurement of the remaining theoretical parameters. 

4.2 Are the NP strong phases negligible? 

The time-dependent angular analysis also allows us to test the assumption that the 
NP strong phases are negligible. Assume that the B — > V1V2 amplitudes contain 
(helicity-dependent) NP strong phases A q x . In this case Eq. (31) can be written 



4 A = Amp(B - 


+ v 1 v 2 ) x -- 


= a x e i5x + A q x e l ^e iA * , 


A x = Amp(B - 


-ViV^ = 


= a A e^+^ A e"^e lA ' 



(44) 

There are now 13 theoretical parameters, and only 11 observables. However, since 
the expressions for the observables in terms of parameters are nonlinear, one can 
obtain bounds on the various theoretical quantities [15]. In particular, one can 
constrain the two NP strong phase differences A^ — Af, % — 0, ||. If either of these 
bounds is inconsistent with zero, this will invalidate the assumption of negligible NP 
strong phases. 



5 Discussion 

In Ref. [8] it was argued that the NP strong phases are negligible relative to those 
of the SM. In this case, all NP matrix elements for a given b — > sqq process (q = 
u, d, s, c) can be summed into a single effective NP operator, with amplitude A q 
and corresponding weak phase $ g . In the previous sections, we have examined three 
different methods for measuring these NP parameters. They are: (i) the comparison 
of b — > s and b — » d penguin B decays [8] , (ii) the combined measurement of B — > 7rK 
and B — > tttt decays, and (iii) the time-dependent angular analysis of B — > V\ V2 
decays. All three methods have their particular advantages and disadvantages. 

There are several B decay pairs to which the -B-penguin method can be applied. 
Depending on which pair is used, the s-, d- and c-quark NP parameters can be 
obtained with a theoretical error of about 5-15%. The w-quark NP parameters A u 
and $ u can also be measured, but with a much larger theoretical error. A key 
assumption is that all SM weak phases have already been measured. Also, most 
pairs of modes involve time-dependent B° s decays, which are hard to measure. 

The B — > tiK J B — >■ tttt method probes only the d- and w-quark NP parameters, 
and with a large (> 25%) theoretical error. Still, the theoretical error on A u and 
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$ u might well be smaller than in the method with B penguin decays. The main 
advantage of this method is that the analysis can be done now. Ref. [31] analyzed 
these B decays within the SM; it is straightforward to include NP parameters in the 
analysis. It is also possible to adapt the method to include NP effects in B — ► titt if 
one makes some assumptions about the nature of the NP. 

The B — > V1V2 method has no theoretical error - flavour SU(3) is not required, 
and we make no assumption about new physics in b — > d transitions. However, this 
method only allows us to measure the s-, d- and c-quark NP parameters, and it 
applies only to those new-physics models in which the NP weak phase is universal. 
Furthermore, it is very difficult experimentally. On the other hand, it can be used 
to test the assumption of negligible NP strong phases. 

At several points previously, we have argued that it is important to measure all 
the NP parameters, and in as many different ways a possible. In the following, we 
show how different NP models lead to different patterns of NP parameters. Thus, the 
measurement of the NP parameters can rule out certain models and point towards 
others. 

There has already been a great deal of theoretical work discussing various NP 
models which can explain the apparent discrepancy in the Belle measurement of 
sin 2/3 in B®(t) — > <f)K s [34, 35, 36]. Our aim here is not to produce an exhaustive 
analysis of such NP models. Instead we consider only two, and show that the 
measurement of the NP parameters can distinguish between them. 

5.1 Z-mediated FCNC's 

The first new-physics model we consider is zT-mediated (or zT'-mediated) FCNC's 
[10]. This model has received much attention as a potential explanation of the 
Bj(t) — > 4>K S result [34]. The Zbs FCNC coupling which leads to the b — »■ s 
transitions is parametrized by the independent parameter Uf b : 

£ Z FCN c = -TT^IT U * trfM* ■ ( 45 ) 

Z COS u w 

Note that the FCNC involves only left-handed s and b quarks. These couplings are 
effectively new contributions to the electroweak penguin operators of the SM. 

The new-physics weak phase arises because U^ b can be complex. However, 
because this parameter is universal, the weak phase of all NP operators will be 
the same. This model therefore predicts the equality of all NP weak phases $ g 
(q = u, d, s, c). If this condition is not found to be satisfied, the model is ruled out. 
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5.2 Supersymmetry with .R-parity breaking 

In supersymmetric (SUSY) models with i?-parity breaking, the relevant part of the 
i?-parity breaking piece is given by 



1 

2' 



Vlfc = K jk LiQiD c k + -X^U'DpZ . (46) 



Here Li and Qi are the left-handed lepton and quark doublet superfields, respectively, 
and Ui and Di are the left-handed quark singlet chiral superfields, where i, j, k are 
generation indices and c denotes a charge conjugate field. 

In the above, the A' and A" couplings violate lepton number and baryon number, 
respectively. The non-observation of proton decay imposes very stringent conditions 
on the simultaneous presence of both couplings [37]. One therefore assumes the 
existence of either L-violating couplings or 5-violating couplings, but not both. 

Now, A"r- fc i is antisymmetric in the last two indices. Thus, there are no B- 
violating couplings which can lead to the b — > sss decay necessary to explain the 
B®(t) — > <pK s result of Belle. We will therefore concentrate only on L- violating 
couplings. (For a discussion of i?-parity violation and B^(t) — *■ 4>K S , see Ref. [35].) 
In terms of four-component Dirac spinors, the L-violating couplings are given by 
[38] 

c xl = -A^[z>i44 + 44^ + (4)*(^lM 

-e!4< - uid k R ei - (4)*(ei) V L ] + h.c. (47) 

From this Lagrangian, we see that there are i?-parity-violating contributions to 
all b — ► s transitions [39]. There is a single contribution to each of the decays 

b — *• suu and b — > sec: 

A' A'* 
Leff = l ^T^u a -i^ L UpS^^ R b a , 

A' A'* 
L eff = TT^Wflififi S(,w R b a . (48) 



Imi 



For b — > sdd, there are four terms: 

-a 



K ff = ^^d lL d7s lR b+^^d lR d7s lL b 



V A'* A' A'* 
_ »i2 ^ dofltfLdp s^^ R b a * 31 | 21 dal^Rdp sp^ L b a . (49) 

Finally, the relevant Lagrangian for the b — ► sss transition is 

L s = *^s lR s s lL b + ^k - slLS - slRh . (50 ) 
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From the above expressions we can deduce the following predictions of ^-parity- 
violating SUSY models. First, in general, all four NP parameters are present and are 
unrelated to one another. Second, since there is only a single term contributing to 
each of b — ► suu and b — > sec transitions, the measured values of $ u and $ c should 
be independent of the decay pairs considered. On the other hand, since there is 
more than one contribution to both b — » sdd and b — > sss, the values of $d and $ s 
will in general be pro cess- dependent. Also, these will differ from $ u and $ c . Should 
this pattern of NP weak phases not be found experimentally, we can either rule out 
or constrain this model of new physics. 

The above two examples illustrate that, indeed, different NP models lead to 
different patterns of the NP parameters A q and $ g . Thus, the knowledge of the NP 
parameters will allow us to discriminate among various models, and rule certain ones 
out entirely. In order to (partially) identify the NP, it will be important to measure 
its parameters in as many different ways and decay modes as possible. A complete 
analysis of NP parameters would therefore include their measurement using all three 
of the methods described in the previous sections. 

6 Conclusions 

The main purpose of the study of CP violation in the B system is to look for physics 
beyond the SM. There are now many theoretical signals of such new physics, and in 
fact there are several experimental hints of NP in decays involving b — ► s penguin 
amplitudes. Still, the conventional thinking was that the identification of NP could 
only be done at future high-energy colliders, where the new particles could be directly 
produced. However, recently it was shown that it is possible to measure the NP 
parameters in B decays [8] . The key observation is that the strong phases associated 
with the NP operators are negligible compared to those of the SM. In this case, all 
NP matrix elements for a given b — ► sqq process (q = u, d, s, c) can be summed into 
a single effective NP operator, with amplitude A q and corresponding weak phase 
$ ? . These NP parameters can be measured, allowing a partial identification of the 
new physics. 

In this paper we have discussed three methods of measuring the NP parameters. 
In most cases, it is assumed that the new physics contributes only to decays with 
large b — > s penguin amplitudes, while decays involving b — > d penguins are not 
affected. The first method, initially proposed in Ref. [8], employs the comparison 
of time-dependent b — » s and b — *■ d penguin B decays. The second uses the 
combined measurement of B — ► nK and B — ► tttt decays. The third requires the 
time-dependent angular analysis of B — > V{Vi decays. 

The three methods can be used to probe different NP parameters, and with 
different theoretical errors. The 5-penguin method allows us to obtain the s-, d- 
and c-quark NP parameters with a theoretical error of about 5-10%. The -u-quark 
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NP parameters can also be measured, but with a much larger theoretical error. The 
B — > TtK jB — > 7T7T method probes only the d- and w-quark NP parameters, but with 
a large (> 25%) theoretical error. The B — ► V1V2 method has no theoretical error, 
but it applies only to NP models with a universal weak phase, and it only allows 
the measurement of the s-, d- and c-quark NP parameters. 

The three methods also have different levels of experimental difficulty. Most 
pairs of modes in the 5-penguin method involve time-dependent B® decays, which 
are hard to measure, and the time-dependent angular analysis of B — ► V1V2 decays 
is very difficult experimentally. On the other hand, the B — ► ttK/B — ► tttt method 
can be performed with present data. 

Ideally, a full analysis of NP parameters would use all three methods. Then all 
the NP parameters A q and § q can be measured in many different ways, and using 
various decay modes. In general, different NP models lead to different patterns 
of the NP operators. The knowledge of the NP parameters will thus allow us to 
discriminate among various models and partially identify the new physics, before 
the direct production of new particles at high-energy colliders. 

Acknowledgements: The work of A.D., M.I., D.L. and V.P. was financially sup- 
ported by NSERC of Canada. 
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u fl - 


= Qal^Lf^f^Lba , 


0% = q^Lff^Lb , 


O q 


= q lfl Lbq'rL(R)q' , 


0% = q alll Lbpq'^L(R)q' a , 


n q 


= ^q llM Lbe q ,qYR(L)q' 


3 

> °8,io = ^l^Lbpeq.q'^RiL)^ , 



Appendix 

In this Appendix, we present a more detailed discussion of the claim that the NP 
strong phases are negligible. To understand our argument, it is crucial to examine 
where the SM and NP strong phases come from. Note that this argument is not 
dependent on any particular calculational framework. However, it will be useful to 
demonstrate it within a particular approach for calculating nonleptonic decays. We 
will follow the method of QCD factorization [25], also known as the BBNS approach. 
To be more specific let us condider B — ► Kn decays. The starting point for 
calculations of this decay is the SM effective hamiltonian for B decays [40]: 

(~i 10 

H q eff = -^[v fb v; q ( Cl oi f + c 2 oi f ) - J2(v tb v t * qC t)o?) + h.c. , (a. i) 

where the superscript t indicates the internal top quark, / can be the u or c quark, 
and q can be either a d or s quark. The operators Of are defined as 



(A. 2) 



in which R(L) = 1 ± 75, and q' is summed over u, d, s, c. O2 and 0\ are the tree- 
level and QCD-corrected operators, respectively. 3 _ 6 are the strong gluon-induced 
penguin operators, and operators O 7 _i are due to 7 and Z exchange (electroweak 
penguins), and "box" diagrams at loop level. 

First, we consider the SM alone. In particular, we consider the QCD penguin 
operators O^q which contribute dominantly to the SM penguin amplitude for B — > 
Kn. (Note that V t bV* s in H e ff can be expressed in terms of V* b V cs and V* b V us 
using CKM unitarity.) Consider the O^q operators with weak phase V* b V cs . (The 
same arguments apply to contributions with weak phase V* b V us .) The SM penguin 
amplitude can obtain a strong phase from two types of rescattering, which we call 
class I and class II: 

• In class I, there is rescattering from the tree-level operators with Wilson coef- 
ficients C\ ~ 1. The operator has the same weak phase as O^q and leads to a 
b — > ccs transition. This will contribute to a final state like Kn only through 
rescattering and will consequently generate a strong phase. The matrix ele- 
ment in question is (Ktt\Oi\B) , and is called P c . In the BBNS picture this 
rescattering is represented by the penguin function G(x) and is proportional 
to C\ ~ 1. The G(x) generates the dominant strong phase. 

In the usual approach, including that of BBNS, the contribution of the tree 
operator is combined with the matrix element of the O^q operator in the 
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quantity 04,6 to represent the SM penguin contribution. The scattering from 
the tree operator is the dominant rescattering effect and is responsible for the 
strong phase of the SM penguin amplitude. Note that the strength of the 04,6 
operators, C^q, is only ~ 4%. Thus, even though rescattering costs a s , it is 
enhanced by the large Wilson coefficient C±, which is 25 times as large as C4. 
The rescattered tree amplitude is therefore roughly the same size as the matrix 
element of the O^q operator, so that it can generate a significant strong phase 
in the BBNS approach [25]. 

• In class II, one can have rescattering from the operators O^q themselves: 
(Ktt\0^q\B). This is represented by BBNS's g(x) and G(x), and is propor- 
tional to C4 or Cq. However, as mentioned above, the size of C4 and Cq is only 
about 4%. Unlike Class I rescattering, there are no operators with large Wil- 
son coefficients to produce an appreciable strong phase here. In other words, 
the strength of these operators is much smaller than C\ ~ 1. So this class of 
rescattering is subdominant. 

Now we consider the presence of NP. At mj,, one can separate the effective 
Hamiltonian into two pieces: O4 (say), which is the SM operator with the SM 
Wilson coefficient, and 4 which contains all the NP contributions. We assume 
that O4 and O4 are of similar size. The strong phase now comes from three sources: 
rescattering from 0\ (class I), O4 (class II) and O4 (class II). The rescattering effects 
from 0\ and O4 are simply those corresponding to the SM, whose strong phase can 
be large. However, the strong phase of the matrix element of O4 comes only from 
rescattering from the NP operator itself: (Ki:\Oi\B). That is, the NP rescattering 
is only of the Class II type. Hence, the NP strong phases are subdominant. To a 
first approximation, we therefore ignore the NP strong phases compared to the SM 
strong phases. 

To demonstrate this point numerically let us work with a specific model of new 
physics, R-parity violating SUSY, that was considered in Sec. 5. In this model the 
Lagrangian for b — > suu and b — *■ sec transitions is : 

A' A'* 
A' A'* 

Kff = T^T^Cal^lLCp S^^ R b a . (A. 3) 

For the decays B — > Kir, the term in L u e y * can contribute directly or through rescat- 
tering of the -u-quark loop in the BBNS method. The quantity r, which is the ratio 
of the imaginary part of the amplitude relative to its real part, is given by [25] 

C F a s {m h ) » 
T = AnN c /m ^ s ")] 
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f 1 
G K {s) = / dxG(s u ,l — x)$£(x) 

Jo 



G(s,x) = -4/ u(l - u)ln[s - u(l - u)x] , (A. 4) 

Jo 

where Cf = (N% — 1)/(2N C ), s u = (m u /mb) 2 and $5(a;) is the light-cone distribu- 
tion(LCD) of the kaon . Using the asymptotic LCD and the fact that s u = to a 
very good approximation, we obtain 

r = \ aK bJ — ^0,015 (A. 5) 

4vrA c 3 v ; 

where we have taken a s (mb) = 0.2. Clearly the NP strong phase from Eq. (A. 5) is 
~ 1°, and is negligible. The factor of a s (mb) / (4:Tt) ps 0.016 comes from the strong 
coupling constant and the loop factor, and produces a significant suppression of the 
rescattering. Only in the SM can this rescattering be significant, as it is propor- 
tional to the large Wilson coefficient C\. There are also rescattering effects such as 
the vertex contributions in the BBNS approach, which will generate new operators 
with different colour structure than !/"//■ However we expect these operators to 
be suppressed by the factor of a s {mb)/{^) ps 0.016 relative to L^jj- and hence to 
be negligible. The Lagrangian L^rr can contribute to Kit final state only through 
rescattering and will be suppressed by the same factor a s {mb)/{^Ti) relative to £"//■ 
It can therefore also be neglected. 

For b — ► sdd, there are four terms in the Lagrangian: 

L d eff = ^^d lL ds lR b+^^d lR ds lL b 
V V* V V* 

%\2 il3 7 t — 7 7'S1 ?21 T 7 — 7 / \ n\ 

- 2 2 daj^Ldp s/WtfRba 2 d^^Rdp s^^ L b a . (A. 6) 

Vi Vi 

For the decays B — > Kn, the terms in L d e ft can contribute directly or through 
rescattering of the ci-quark loop in the BBNS method. The terms involving vec- 
tor operators will have the same rescatterieng as L u e ^ and so the resulting strong 
phases can once again be neglected. The scalar operators can only have electroweak 
rescattering though the the exchange of a photon or a Z and will therefore be even 
smaller. 

Finally, the relevant Lagrangian for the b — > sss transition is 

y v* \' \'* 

ts A i32 A i22 - - 1 . A i22 A i23 - - 1 I A 1\ 

L eff = 9 S ^ S S 7Lb H n S ^ LS S 7R & • ( A - <) 

11 m~ mA 

Vi v t 

The terms in L s e jr can contribute to the K7r final state only through rescattering. 
They will therefore be suppressed by the factor a s (m&)/(47r) relative to £"//' and 
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hence can be neglected. Note that the operators sbbb and shit can also contribute 
through rescattering but will have no imaginary contribution as the gluon momen- 
tum is below the threshold of bb and tt production. 

We have therefore shown in general terms, and within a specific model of NP, 
that the NP strong phases are small compared to those of the SM, using the BBNS 
approach. As mentioned above, this result is independent of the specific model of 
NP, as well as of the method used to calculate nonleptonic decays. 
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